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Zusammenfassung 



1 Bezeichnung 

Vorrichtung und Verfahren zur Bewertung von Finanzderivaten mit Hilfe von diinnen Gittern. 

2.1 Technische Aufgabe und Zielsetzung 

Die schnelle und genaue Bewertung komplexer Finanzderivate ist ein zentrales Problem vieler 
Banken, Versicherungen und Anleger. In vielen Fallen kann der Preis solcher Derivate als Erwar- 
tungswert bestimmt werden. Bestehende Vorrichtungen und Verfahren zur Berechnung dieser 
Erwartungswerte basierend auf Monte Carlo oder Quasi-Monte Carlo Integration sind jedoch 
sehr zeitintensiv und ungenau. 

2.2 Losung der Aufgabe 

Die Bewertung von Finanzderivaten wird mittels einer Vorrichtung bestehend aus einem Com- 
puter und einem Verfahren basierend auf multivariater numerischer Integration mit Hilfe von 
diinnen Gittern durchgefiihrt. Das vorgeschlagene Verfahren kann die auftretenden Erwartungs- 
werte effizienter als bisherige Verfahren berechnen und damit die Preise der Derivate wesentlich 
schneller und genauer ermitteln. 



Die Vorrichtung beinhaltet eine Ein- (1) 
und eine Ausgabeeinheit (6), eine Spei- 
chereinheit (5) sowie eine Programmein- 
heit (2-4). Die Programmeinheit besteht 
aus einem Aufbaumodul (2), einem Dis- 
kretisierungsmodul (3) und einem Integra- 
tionsmodul (4). 

Im Falle nicht glatter Integranden, wie 
sie z.B. bei der Bewertung von Optio- 
nen auftreten, wird das Integrationsgebiet 
durch Transformationen in glatte Teilge- 
biete zerlegt. Eine zusatzliche Beschleu- 
nigung des Verfahrens wird durch Di- 
mensionsreduktion sowie Parallelisierung 
ermoglicht. 



2 • 3 Anwendungsgebiet 

Die Vorrichtung ermoglicht zusammen mit dem vorgeschlagenen Verfahren die schnelle und 
genaue Bestimmung objektiver sowie subjektiver Preise von Finanzderivaten, wie z.B. Bonds, 
Swaps, Futures, CMOs und Optionen. Die ermittelten Preise konnen z.B. zum Kauf oder Verkauf 
der Derivate selbst, zur Portfoliooptimierung, Preisbestimmung von Vertragen, Risikomanage- 
ment oder Entscheidungsfindung bei Investitionen verwendet werden. 
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Beschreibung 

1 Technisches Gebiet 

Die Erfindung liegt auf dem Gebiet der elektronischen Datenverarbeitung, speziell auf dem 
Gebiet der Finanzinformatik. Insbesondere betrifft die Erfindung die schnelle und genaue Be- 
wertung von Finanzderivaten mit Hilfe von elektronischen Rechnensystemen. 

2 Stand der Technik 

Die Bewertung von Finanzderivaten hat im letzten Jahrzehnt einen starken Aufschwung erfah- 
ren und ist zu einem wichtigen Instrument in der Finanzwirtschaft geworden. Einerseits ist die 
Moglichkeit einer Bewertung eine notwendige Voraussetzung, damit ein Finanzinstitut iiber- 
haupt ein neues Derivat anbieten kann, andererseits ist die Bewertung notig um ein Portfolio 
geeignet zu strukturieren, abzusichern und zu optimieren. * ^ 

Es gibt nun eine Reihe verschiedener Arten von Finanzderivaten, die beispielsweise Zins- * 
oder Aktien-basiert sind. Hierzu zahlen u.a. Bonds, Swaps, Futures, CMOs und Optionen. Ihre 
Bewertung erfolgt, unter der Annahme, dafi keine Arbitragemoglichkeit besteht, mittels parti- 
eller Different ialgleichungen oder iiber den Martingalansatz. Letzterer ist jedoch allgemeiner. 
Dabei spezifiziert man zunachst einen stochastischen Prozefi fur den zugrunde liegenden Wert 
(zum Beispiel den Zinssatz oder den Aktienkurs). Dann bestimmt man das aquivalente Mar- 
tingalmafi, das den zugrunde liegenden Prozefi in ein Martingal umwandelt und berechnet den 
Wert des Derivats als den Erwartungswert seiner discounteten Auszahlungsfunktion unter die- 
sem risikoneutralen Mafi. Unter bestimmten Voraussetzungen ist der Zugang iiber eine partielle 
Different ialgleichung und das Martingal gleichwertig. Der Zusammenhang wird dabei iiber ver- 
allgemeinerte Feynman-Kacz Formeln hergestellt. Der Martingalansatz ist jedoch universeller 
einsetzbar und leichter auf neue Situationen iibertragbar. 

Der Preis eines Finanzderivats lafit sich somit als Erwartungswert ausdriicken. In konti- 
nuierlicher Zeit beinhaltet der Integrand dabei wiederum ein Pfadintegral, das etwa mit dem 
Euler-Verfahren/Trapezregel diskretisiert werden kann oder es wird bedingt durch die Anwen- 
dung gleich ein zeitdiskretes Modell verwendet. In beiden Fallen fuhrt dies zu einem hochdi- 
mensionalen Integrationsproblem. Fur die einfache Europaische Call Option ist es moglich eine^p; 
geschlossene Losung anzugeben (Black-Scholes Formel), kompliziertere Optionstypen benotigen 
jedoch ein numerisches Losungsverfahren. Dies gilt analog auch fur die anderen Arten von Fi- 
nanzderivaten. 

Als numerisches Integrationsverfahren fur hochdimensionale Integranden ist die klassische 
multivariate Quadratur nicht geeignet. Hier begegnet man dem Fluch der Dimensionalitat: Der 
Aufwand skaliert hier exponentiell mit der Dimension. Man hat die Komplexitat der Ordnung 
0(N~ r/d ) wobei r die Glattheit des Integranden und d seine Dimension bezeichnen. Hingegen 
ist das Monte Carlo Verfahren von der Dimension unabhangig. Dabei wird der Integrand / an 
einer zufalligen Folge von N Punkten ausgewertet und man erhalt die Quadraturformel 

1 N 

Das Monte Carlo Verfahren konvergiert jedoch nur sehr langsam (und auch nur im statisti- 
schen Sinne). So ist die Genauigkeit, die mit N Funktionsauswertungen erzielt werden kann, 
von der Ordnung 0(l/y/N). Von besonderem Interesse sind deswegen die im letzten Jahrzehnt 
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entwickelten Quasi Monte Carlo Methoden [4]. Hier wird der Integrand fur eine deterministisch 
bestimmte Folge von Punkten x; ausgewertet und man verwendet als Quadaturformel analog 
zum Monte Carlo- Verfahren 

2=1 

Es gibt hierbei eine Reihe verschiedener Konstruktionen, beispielsweise die Halton-, Sobol- oder 
Faure-Folgen, die sich in ihrer Praasymptotik unterscheiden aber alle in einer Konvergenzord- 
nung des Typs 0((logN) d /N) resultieren. Diese ist also fast eine Ordnung besser als fur die 
Monte Carlo Methode. Zudem ist der Fehler deterministisch. Prototyp ist hierbei das Programm 
FinDer von J. Traub [9], das mittlerweile in vielen Banken verwendet wird. Diese Methode wurde 
in den Patenten US005940810A [7] und US006058377A [8] dokumentiert. 

Weiterhin gibt es sogenannte Dunngitterverfahren [2, 6]. Bei diesem Ansatz werden mehrdi- 
mensionale Quadraturformeln durch geeignete Kombination von Tensorprodukten eindimensio- 
naler Quadraturformeln, wie z.B. Clenshaw-Curtis oder Gaufi-Patterson Formeln, konstruiert. 
Das allgemeine Dunngitterverfahren kann wie folgt beschrieben werden. Man betrachtet zunachst 
eine Folge eindimensionaler Quadraturformeln fur eine eindimensionale Funktion /, 

n] 

Qlf ■= Y. W li -f( X li)- 
i=l 

Nun definiere die Differenzformel durch 

Aj/ := (Ql~Ql-i)f with 
Qlf := 0. 

Die Diinngitterkonstruktion fur d-dimensionale Funktionen / besteht fur I G IN und kEN d aus 

Qff:= ^(Ai^.-.OAJ,)/. 

mit Indexmengen X{ sodafi fiir alle k € Ii, 

k - Gj e Xi for 1 < j < d, kj > 1, 

gilt. Spezialfalle des Verfahrens sind klassische dunne Gitter, wo I« = {|k|i < / + d - 1 }, sowie 
klassische Produktformeln, bei denen T { = {\k\oo < I}- Abb. 4 gibt Beispiele verschiedener klas- 
sischer diinner Gitter im 2D Fall basierend auf der Trapez-, der Clenshaw-Curtis Formel, der 
Gaufi-Patterson und der Gaufi-Legendre Formel. Ein substanzieller Unterschied zu Monte Carlo 
und Quasi-Monte Carlo Verfahren ist hier die Verwendung unterschiedlich grofier Gewichtsfakto- 
ren, d.h. in der Darstellung Q N f = Eili mf(*i) sind die Gewichte Wi des Dunngitter verfahrens 

nicht alle gleich l/N. (d-D(r+i) 

Die Konvergenzordnung das klassischen Dunngitterverfahrens ist e = 0{\og{NY ^ r >N r ) 
und damit ebenfalls unabhangig von der Dimension d. Im Gegensatz zu Monte Carlo und Quasi- 
Monte Carlo Verfahren kann das Dunngitterverfahren allerdings auch die Glattheit r des Inte- 
granden nutzen und hat damit exponentielle Konvergenz bei glattem Integranden (r -> oo). Fiir 
glatte Integranden ist damit dieses Verfahren substantiell schneller als die Monte Carlo oder 
Quasi-Monte Carlo Verfahren. 

Fur das Dunngitterverfahren ergeben sich jedoch zwei grundsatzliche Probleme: Zum emen 
verschlechtert sich die Konvergenzrate des Verfahrens bei nicht glatten Integranden, wie sie 
z.B. bei Optionen auftauchen, d.h. es verliert in diesem Fall seinen Vorteil. Zum anderen ist 
das Verfahren ebenso wie Quasi-Monte Carlo Verfahren nicht vollkommen unabhangig von der 
Dimension und die Konvergenzrate verschlechtert sich ebenfalls mit ansteigender Dimension. 



3 Erfindung 



Die Erfindung bezieht sich hierbei auf eine Vorrichtung bestehend aus einem Computer und ein 
Verfahren, welches auf der Dunngittertechnik basiert, die die schnelle Berechnung von Erwar- 
tungswerten erlaubt, wie sie zur Bewertung von Derivaten aus dem Finanz- und Bankenwesen, 
z.B. Zins-, Aktien-, Wahrungs- oder Giiterderivaten, benotigt werden. 

Die Vorrichtung (Abb. 1) besteht aus einem Computer mit einer Ein- (1) und Ausgabeeinheit 
(6), einer Speichereinheit (5) sowie einer Programmeinheit (2-4). Die Ein- und Ausgabeeinhei- 
ten ermoglichen die Eingabe der Derivatparameter und die Ausgabe des Derivatwertes. Die 
Speichereinheit wird zur Ablegung des Programms, der Ein- und Ausgabewerte, sowie von Zwi- 
schenergebnissen benotigt. Die Programmeinheit besteht aus einem Aufbaumodul (2), einem 
Diskretisierungsmodul (3) und einem Integrationsmodul (4). Das Aufbaumodul hat die Aufgabe 
der Bestimmung der Integrationsfunktion. Das Diskretisierungsmodul ermittelt geeignete Inte- 
grationspunkte mittels diinner Gitter. Das Integrationsmodul berechnet daraufhin den Wert des 
Derivats durch Kombination der Funktionswerte an den Integrationspunkten. 

Im Folgenden werden die einzelnen Komponenten genauer beschrieben. Die Eingabe erfolgW 
liber eine Tastatur, ein anderes Eingabegerat oder iiber eine analoge bzw. digitale Verbindung. * ' 
Die Eingabedaten spezifizieren den genauen Typ des Derivats sowie die Eigenschaften der zu- 
grunde liegenden Werte. Die Daten sowie weitere Zwischenergebnisse werden in einem Speicher 
abgelegt. Ein Computerprogramm liest in einem ersten Schritt (Aufbau) diese Daten ein und 
baut ein mehrdimensionales Integrationsproblem basierend auf den Eingabedaten auf. In einem 
zweiten Schritt (Diskretisierung) werden mit Hilfe dieser Informationen geeignete Integrations- 
punkte und Integrationsgewichte basierend auf der Dunngitter-Technik ermittelt. Im dritten 
Schritt (Integration) wird das Integrationsproblem durch Auswertung des Integranden an den 
berechneten Integrationspunkten und Kombination der Ergebnisse gelost und als Ergebnis der 
Wert des Derivats zuriickgeliefert. Die Ausgabe erfolgt iiber einen Bildschirm, ein anderes Aus- 
gabegerat oder eine analoge bzw. digitale Verbindung. 

Dieses Verfahren ist bei glatten Integranden, wie sie typischerweise zur Bewertung von Zins- 
derivaten auftreten, ohne weitere Modifikationen einsetzbar. Beispiele hierfur sind Portfolio- 
Optimierung, Preisfestsetzung von futures und forward- Vertragen, die Bewertung hypotheken- 
basierter Wertpapiere, Preisbestimmung von Versicherungsvertragen, Risikomanagement oder 
Entscheidungsfindung bei Investitionen. Im Falle von nicht-glatten Integranden, wie sie z.B. bei 
der Bewertung von Optionen auftreten, verschlechtert sich jedoch Konvergenz des Verfahren?^ , 
stark. In so gut wie alien Fallen besitzen hier die Integranden unstetige erste Ableitungen (r = 1), 
in manchen Fallen ist sogar der Integrand selbst unstetig (r = 0). Ein wesentlicher weiterer Be- 
standteil der Erfindung ist daher eine Zerlegung des Integrationsgebiets in glatte Bereiche im 
Diskretisierungsschritt. Durch geeignete Transformationen werden dann die Diinngitterquadra- 
turformeln auf diese Teilbereiche abgebildet und das Gesamtintegral als Summe dieser Einzelin- 
tegrale berechnet. Auf diese Weise wird nur iiber glatte Bereiche integriert und die positiven 
Eigenschaften des Dunngitterverfahrens bleiben erhalten. 

Das Dunngitterverfahren ist zwar weitestgehend, aber dennoch nicht vollstandig von der Di- 
mension des Integrationsproblem unabhangig. Die Dimension d tritt in der Konvergenzordnung 
als Exponent eines logarithmischen Faktors auf. Dies fuhrt zu einer (wenngleich auch langsamen) 
Verschlechterung der Konvergenz, wenn die Dimension d ansteigt. Daher ist es sinnvoll, dimensi- 
onsreduzierende Verfahren anzuwenden. Zum einem wird im Falle von pfadabhangigen Optionen 
der zugrunde liegende stochastische Prozesses hierarchisch (z.B. mit Braunscher Briicke) diskre- 
tisiert. Zum anderen wird im Falle von performance-abhangigen Optionen die Volatilitatsmatrix 
mittels Singularwertzerlegung transfer miert. In beiden Fallen wird durch eine Fokussierung auf 
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die jeweils wichtigsten Dimensionen im Aufbaumodul eine Reduktion der effektiven Dimensi- 
on vorgenommen. Das Dunngitterverfahren verwendet dann im Diskretisierungsmodul genauere 
Quadraturformeln in wichtigeren Dimensionen und weniger genaue Formeln in weniger wichti- 
gen Dimensionen. Auf diese Weise wird in vielen Fallen eine weitestgehende Unabhangigkeit des 
Verfahrens von der Dimension und in den restlichen Fallen zumindest eine starke Beschleunigung 
des Verfahrens erreicht. 

Eine zusatzliche Beschleunigung des Verfahrens geschieht durch Parallelisierung, d.h. ei- 
ne verteilte Berechnung auf mehreren Prozessoren. Sie ist ein weiterer Bestandteil der Erfin- 
dung. Hierbei wird das gesamte Integrationsproblem in voneinander unabhangige Teilprobleme 
unterteilt. Dies sind z.B. die einzelnen Teilbereiche, die durch obige Transformation entstan- 
den sind, oder die Teilsummen iiber die verschiedenen Multiindizes k aus der Diinngitter- 
Quadraturformel. Die Vorrichtung ist in Abb. 2 und 3 dargestellt. Aufgabe des Verteilungs- 
moduls (7) ist die Zerlegung des Integrationsproblems. Die Kombination der Teilergebnisse wird 
wahrend der parallelen Rechnung durch ein Zusammenfuhrungsmodul (8) erledigt. Die einzel- 
nen Teilprobleme werden entweder in einem parallelen Rechensystem mit verteiltem Speicher 
(Abb. 2) oder gemeinsamen Speicher (Abb. 3) durchgefuhrt. 

4 Ausfuhrungsbeispiel 

Eine Option ist ein Vertrag zwischen zwei Parteien, der dem Halter das Recht aber nicht die 
Pflicht gibt, eine festgesetzte Menge von z.B. Aktien zu einem bestimmten Zeitpunkt zu einem 
vorher festgelegten Preis zu kaufen/verkaufen. Eine Option stellt damit einen intrinsischen Wert 
dar, da der Halter die Option nicht ausiiben mufi. Die Frage ist nun wie der Preis fur eine Option 
fair bewertet werden kann. Konkret haben wir hier mit dem Martingalansatz die Bewertung 

V = e-r T E(P({S{})). 

Dabei bezeichnet r die konstante Zinsrate, t = T den Zeitpunkt der Optionsausubung, ]P die 
Auszahlungsfunktion, S{ den Aktienwert der i-ten Aktie des Portfolios zur Zeit jAt, mit At = 
T/M. Hier ist M die Zahl der diskreten Zeitpunkte j = 1, „,M und N bezeichnet die Zahl der 
betrachteten Aktien % = l,..,iV. e~ rT ist der Discountfaktor und t = 0 ist der Zeitpunkt der 
O p t ionsb e wer t ung . 

Wir nehmen nun an, dafi die Aktien des Portfolios einem Modell folgen, das durch die 
stochastische Differentialgleichung 

N 

dSi = frSidt + ^2 °ikS k dW k (1) 

beschrieben ist, wobei (a ik ) die zeitkonstante Volatilitatsmatrix bezeichnet, ^ die zeitkonstante 
Drift fur die Aktie % bezeichnet 1 und dWi ein geometrischer Brownscher Prozefi fiir die Aktie i 
ist, d.h. Wi ist das Wiener Mafi. Integration und Ito-Formel ergibt nach Diskretisierung in der 
Zeit mit den Zeitpunkten jSt die Werte Si 

5 j+i = S 3 . e (^-|(ELi-? fc )) At +Er=i^^^' (2) 

und somit gilt 2 

S i = 5 o e Ef =1 [(^-i(ELi^)) At +ELi^^^]. 

1 Es gibt auch Ansatze mit zeitabhangigen Volatilitatsmatrizen und zeitabhangiger Drift. Dann folgen a ik und 
fj,i eigenen SDEs, die ihr Verhalten modellieren. 

2 Im Computerprogramm wird die erste rekursive Formel genutzt. Die jetzt folgende nicht-rekursive aufsum- 
mierte Form ist aber notwendig um den spateren Integranden definieren zu konnen. 
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Hier sind W l k Zufallsvariable, die iV(0 ; l)-verteilt sind. 

Die Definition des Erwaxtungswerts E(.) ist nun fur eine allgemeine Funktion / 

/oo roo 
•■•/ Mx*}) g{x\)...g{x%) dx\...dx 
-oo J — oo 

mit der Standard-Normalverteilung 

- 

Mit der Forderung, dafi keine Arbitrage existiert, mufi der Erwartungswert beziiglich des zum 
unterliegenden Prozesses aquivalenten Martingalmafies genommen werden. Es werden dabei in 
unserem Fall im Integral des Erwartungswerts die m durch r ersetzt. Damit findet auch der 
Wechsel der SDE (1) in die risikoneutrale Form statt. 

Was nun noch fehlt sind die konkreten Auszahlungsfunktionen. Diese sind abhangig von 
jeweiligen Typ der Option. Generell haben sie die Struktur 

P({S{}) = max(0, H({S{}) - K) (3)^ 

fur Call Optionen und 

P({Si}) = uua{O i K-H({S{})) 

fur Put Optionen. Dabei ist K der Ausiibungspreis. Beispiele sind pfadabhangige Optionen 
(N = 1, M > 1) und performanceabhangige Optionen (N > 1, M — 1). 

Popularstes Beispiel fur pfadabhangige Optionen sind sogenannte Asiatische Optionen. Hier- 
bei wird fur die Auszahlungsfunktion der Mittelwert iiber die Aktienkurse zu alien Zeitpunkten 
zwischen dem momentanen Zeitpunkt und dem Ausiibungszeitpunkt genommen. Wir betrachten 
im folgenden den Fall des geometrischen Mittels, d.h. H({S{}) = (Sj - S f • . . . • Sf*) 1 /™ . Hierzu 
existiert eine geschlossene Losung in Form einer verallgemeinerte Black-Scholes Formel. Kleine 
Anderungen und Variationen (z.B. schon das arithmetische Mittel) sind jedoch nicht mehr ohne 
weiteres analytisch losbar und benotigen ein spezielles numerisches Integrationsverfahren. 

Ein Beispiel fur performanceabhangige Optionen sind sogenannte "Tailored Options" [3]. 
Hier hangt der Preis der Option von der relativen Performance eines Aktienkurses in bezug auf 
die anderen Aktienkurse einer betrachteten Menge von Aktien ab. Ein Beispiel fur eine Call 
Option ist /> 

prKn , _ / max(Si - K t 0) falls Sj/S? > S}/Sl 1 < i < N 

i)-^ a , max (5i _ K% o) sonst ^ } 

wobei a ein Parameter aus [0, 1) ist. Eine Mischung aus beiden Optionstypen ist selbstverstand- 
lich auch moglich. Mit M = N = 1 ergibt sich als Spezialfall die Europaische Call Option mit 
der Black Scholes Formel als analytischer Losung. 

Nun ist also zur Optionsbewertung ein im allgemeinen hochdimensionales Integrationspro- 
blem zu losen. Die Dimensionen riihren dabei von den Zeitpunkten her, liber die bei pfad- 
abhangigen Call Optionen gemittelt wird, oder/und von der Zahl der betrachteten Aktien bei 
performanceabhangigen Optionen. Aufgrund der Definition des Erwartungswerts ist das Inte- 
grationsgebiet bisher noch ( — 00,00)^"^. Um Integrationsmethoden anwenden zu konnen trans- 
formieren wir das Integral des Erwartungswerts mit Hilfe der kumulativen Normalverteilung 
G(y) := f^ OQ g{x)dx auf [0, 1]^ M und erhalten damit 

E(f({x{})) = /\ . . /' f({G(x{)}) dx\... dx%r. (5) 
JO Jo 
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Fur die Integration von (5) mit (3) bzw. (4) wollen wir nun die Dunngitter-Technik gewinn- 
bringend einzusetzen. 

Fur Optionen ergibt sich nun jedoch folgendes Problem: Die Auszahlungsfunktion ist durch 
die Natur einer Option nicht mehr glatt. Dies driickt letztendlich aus, dafi man eine Option ja 
nicht ausuben wird wenn der Kauf oder Verkauf der abgesicherten Aktie einen Verlust ergeben 
wurde. Der Integrand weist bzgl. einer (M - N - l)-dimensionalen Mannigfaltigkeit einen Knick 
(pfadabhangige Optionen) oder sogar Sprunge (performanceabhangige Optionen) auf 3 . Beispiele 
fur den Integranden im zweidimensionalen Fall sind in Abb. 5 gegeben. Deutlich sieht man den 
Knick bzw. den Sprung des Integranden. Die Glattheitsvorraussetzungen fur die Anwendung des 
Dunn-Gitter-Verfahrens sind damit nicht mehr gegeben. 

Die Idee ist nun nur iiber den Trager des Integranden zu integrieren. Hier ist der Integrand 
eine glatte Funktion, die Knicke und Sprunge befinden sich gerade am Rand des Tragers. Urn 
dieses Gebiet zu bestimmen geniigt es, die Nullstellen des Integranden zu berechnen. Wenn man 
das Integral iteriert berechnet reduziert sich die Nullstellensuche auf eine einzige (die letzte) 
Dimension. Wir bestimmen deswegen in der letzten Dimension die Nullstelle x (Newtonverfahren 
fur den Knick oder Bisektionsverfahren fur den Sprung) und transformieren den Integranden 
bezuglich der letzten Dimension mit der linearen Abbildung 

t(x) = x • (1 - x) 4- x 

auf [0,1]. 

Abb. 5 gibt einen Vergleich der verschiedenen Integrationsverfahren bei der Bewertung einer 
pfadabhangigen Option mit 6 Zeitpunkten, Af = 6,JV = 1, und einer performanceabhangigen 
Option mit M = l,iV = 2. Die iiberlegene Konvergenzordnung unseres Diinn-Gitter Verfah- 
rens mit Transformation und Gaufi-Patterson Formel iiber die anderen Methoden (Monte Carlo 
(MC), Klassischer Produktansatz ohne (PR) und mit Transformation auf den Trager (PRTR), 
Quasi-Monte Carlo ohne (QM) und mit Transformation auf den Trager (QMTR), sowie Gaufi- 
Patterson-Dunne-Gitter ohne (SG) und mit Transformation auf den Trager (SGTR)) ist deutlich 
zu sehen. Man beachte, dafi der Fehler hier logarithmisch aufgetragen ist. Es zeigt sich, dafi un- 
ser neues Verfahren den Monte Carlo und Quasi-Monte Carlo Methoden iiberlegen ist. Diese 
konnen eine hohere Glattheit des Integranden nicht ausnutzen, die neue allgemeine Diinngitter- 
Quadratur kann dies hingegen in optimaler Weise. M-iHr+D 

Nachdem die Komplexitat des Dunngitterverfahrens aufgrund des Terms log(JV) 1 Mr > 
nicht vollkommen unabhangig von der Dimension des Problems ist, ist es sinnvoll, zusatzlich 
Dimensionsreduktionsverfahren einzusetzen. Hier soil erlautert werden, wie das Verfahren durch 
hierarchische Diskretisierung des stochastischen Prozesses und adaptive Verfeinerung das Ver- 
fahren weiter beschleunigt werden kann. 

Der naturlichste Weg einen stochastischen Prozess zu diskretisieren ist durch einen random 
walk, d.h. durch die rekursive Formel 

wobei 6(W/) genau dem Exponenten aus der Formel 2 entspricht. Durch Diskretisierung mittels 
Braunscher Briicke [1], wird jedoch der Prozefi mittels eines zukunftigen und eines vergangenen 
Wertes diskretisiert 

gj = S i + S i c bU/kAt/2Wj) 

^ 2 

3 Neben den Singulaxitaten am Rand des Integrationsgebiets, die durch die Transformation auf [0, l] d entstanden 
sind. 
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So werden startend mit S t °, S? := 2^^) die Werte Sf ' 2 , S t M/4 , S? M/ \ S t M/8 , S™'*,... 
usw. ermittelt. Dies fuhrt zu einer Konzentration der Gesamtvarianz des Prozesses in den er- 
sten Schritten der Diskretisierung was die Konvergenzrate von Quasi-Monte Carlo Methoden 
verbessert. 

Fur das klassische Dunngitterverfahren ergibt sich kein sofortiger Vorteil durch diese Dis- 
kretisierungstechnik nachdem alle Dimensionen von gleicher Wichtigkeit sind. Das allgemeine 
Dunngitterverfahren kann jedoch nun dimensionsadaptiv eingesetzt werden und Quadraturfor- 
meln mit niedrigerem Grad in weniger wichtigen Dimensionen verwenden. Auf diese Weise re- 
duziert sich die effektive Dimension fur Integrale vom Typ (5) bei pfadabhangigen Derivaten, 
wodurch die entstehenden Integrationsprobleme schneller und genauer berechnet werden konnen. 
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Abb. 6 



Patentanspriiche 



1. Vorrichtung und Verfahren zur Bewertung von Finanzderivaten, wobei der Wert des Deri- 
vats durch die Ermittlung eines Erwartungswerts berechnet wird, bei dem 

(a) mittels einer Tastatur oder eines anderen Eingabegerats die Parameter des Derivats 
eingegeben werden, 

(b) mittels eines Computers basierend auf den Eingabeparametern ein Integrand aufge- 
baut wird und ein geeigneter mehrdimensionaler Integrationsbereich bestimmt wird, 

(c) mittels eines Computers Integrationspunkte und Integrationsgewichte durch ein Dunn- 
gitter- Verfahren ermittelt werden, 

(d) mittels eines Computers der Integrand im Integrationsbereich an den Dunngitter- 
Integrationspunkten ausgewertet wird, 

(e) mittels eines Computers der Erwartungswert durch Kombination der Integrandwerte 
unter Verwendung der Integrationsgewichte berechnet wird. 

(f) mittels eines Bildschirms oder eines anderen Ausgabegerats der berechnete Wert des 
Derivats ausgegeben wird. 

2. Vorrichtung und Verfahren nach (1), wobei Ein- und Ausgabe iiber eine analoge oder 
digitale Verbindung stattfinden. 

3. Vorrichtung und Verfahren nach (1), wobei der Integrationsbereich durch Nullstellensuche 
bestimmt wird. 

4. Vorrichtung und Verfahren nach (1), wobei die Integrationspunkte sowie Integrationsge- 
wichte dynamisch wahrend der Berechnung des Integrals ermittelt werden. 

5. Vorrichtung und Verfahren nach (1), wobei das Integrationsgebiet in mehrere Teilintegrale 
zerlegt wird und mehrere Prozessoren zur Berechnung der Teilintegrale verwendet werden. 

6. Vorrichtung und Verfahren nach (1), wobei mehrere Prozessoren zur Berechnung der In- 
tegrationspunkte, Integrationsgewichte, zur Auswertung des Integranden oder zur Kombi- 
nation der Integrandwerte herangezogen werden. 

7. Vorrichtung und Verfahren nach (1), wobei die Dimension der Integrale zur Beschleunigung 
des Verfahrens reduziert werden. 

8. Vorrichtung und Verfahren nach (1), wobei zur Bewertung mehrere Erwartungswerte be- 
rechnet und kombiniert werden. 



